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We perform density functional theory calculations on a series of armchair and zigzag nanotubes of 
diameters less than lnm using the all-electron Full-Potential(-Linearised)-Augmented-Plane-Wave 
(FPLAPW) method. Emphasis is laid on the effects of curvature, the electron beam orientation 

■ and the inclusion of the core-hole on the carbon electron energy loss K-edge. The electron energy 
loss near-edge spectra of all the studied tubes show strong curvature effects compared to that of flat 
graphene. The curvature induced ty — a hybridisation is shown to have a more drastic effect on the 
electronic properties of zigzag tubes than on those of armchair tubes. We show that the core-hole 
effect must be accounted for in order to correctly reproduce electron energy loss measurements. 
We also find that, the energy loss near edge spectra of these carbon systems are dominantly dipole 
selected and that they can be expressed simply as a proportionality with the local momentum 
projected density of states, thus portraying the weak energy dependence of the transition matrix 

, elements. Compared to graphite, the ELNES of carbon nanotubes show a reduced anisotropy. 

<n : 

m . 

Q , I. INTRODUCTION 

m : ... n . 

Since the discovery of carbon nanotubes by Iiiima [1| in 1991, much effort has been devoted ex- 
perimentally 0, Q, 0, 0] and theoretically 0, Q, Q to study this new material. The one dimensional 
character of single wall nanotubes (SWNT's) enables them to exhibit very interesting physical 

■ properties due to the quantum confinement of electrons in a one dimensional lattice. Remarkable 
optical, thermoelectric, mechanical and electronic properties can be expected. Doping or inter- 
calation Q, E3, EH E3, E3, E3, E3, E3, structural defects Q Q and pressure effects 

, ^ ' are known to enrich these properties. Resonant Raman spectroscopy and optical absorption spec- 



troscopy |2lL 1221 12.1 |24j have been widely used to describe the electronic structure of isolated 
and ropes of nanotubes. Scanning tunneling microscopy has also been used to detect the 

Van Hove spikes of SWNT's. Most of the theoretical studies are based on the application of the 
Born- von Karman boundary conditions to the two dimension al grap hene sheet pi 13 m a zone- 
folding technique and the all- valence tight- binding method 0, hd-h^ . This graphene zone- folding 
approach is known to be inaccurate to describe the unoccupied density of states (DOS) [3 a few 
eV's beyond the Fermi level. Density functional theory (DFT) calculations |2g, |29j have also been 
carried out but very few calculations [3 probing core-loss spectra are available. Electron energy 
loss spectrum (EELS) and X-ray absorption spectroscopy (XAS) measurements are available but 
their theoretical interpretation is often limited to the site and angular momentum projected lo- 
cal density of states (LDOS). This interpretation does not account for transition matrix elements 
which are known to be important and sensitive to the momentum transfer direction for anisotropic 
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materials like graphite and carbon nanotubes. While the DOS of a tube is characterised by the 
Van Hove singularities inherent to the one dimensional nature of the tube, it has been argued that 
the EEL spectrum is free of any such spikes due to the effect of the Cls core- hole in the final 
state 30j. This implies a strong energy dependence of the transition matrix elements. We point 
out, contrary to this assertion, that the energy dependence of the transition matrix elements is 
weak meaning that the spikes may always feature in the EELS spectra if the energy resolution is 
good enough and thermal effects are rendered minimal. The few DFT calculations of energy loss 
near edge spectra (ELNES) so far available J3| do not include the all important core-hole effect 
on the ELNES. We have included and analysed these effects on the carbon K-edges of nanotubes. 

We perform ab initio calculations on a series of single-wall armchair and zigzag nanotubes. We 
present results on the electronic properties of the carbon nanotubes. For comparison with other 
calculations, we considered the (6,6), (5,0) and (10,0) tubes. The obtained band structure is in 



good agreement with earlier results for both the (6,6) tube |28| and the (5,0) tube |31|,|32j. The 
calculated C K-edges agree with those of Suenaga et al.|3| for the (5,0) and (10,0) tubes. This 
agreement with other calculations validate the method used in this work. Much emphasis is laid 
on the curvature, orientation and core-hole effects on the C K-edges, as these aspects have not 
been addressed so far. 

In order to demonstrate the fact that energy loss near edge structure (ELNES) probes the local 
bonding states in materials, we considered so-called nanoarcs, which consist of curved graphite 
sheets which are repeated in space either as wavy objects or as bonded arcs. For feasible calcula- 
tions, we consider the extreme cases of the armchair (n,n) and the zigzag (n,0) nanotubes which 
form chiral angles of 30° and 0° , respectively. All other nanotubes (chiral) have their chiral angles 
within the range [0°-30°]. 

This article is organised as follows. Section 2 is devoted to the description of the theoretical 
background. In section 3 we give the computational details. Section 4 is devoted to a discussion 
of the main results and we end in section 5 with a brief summary. 

II. THEORETICAL FRAMEWORK 
A. The augmented plane-wave method 

Within the density functional theory the ground state properties are derived from the ground 
state density n(r) which minimises the total energy functional 



E tot [n]=T{n] + U[n}+E xc [n], 



(1) 
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where T is the kinetic energy functional of the non-interacting electron system, U is the electrostatic 
energy functional and E xc is that of the exchange and correlation energy. The electron den sity 
which minimises E to t[n] is obtained by solving self-consistently the Kohn-Sham equations |33t l34j 

H^i (r) = [f + V eJ f (')] V'i ( r ) = ( r ) ( 2 ) 

where ej and tpi are the Kohn-Sham eigenvalues and eigenfunctions, respectively. The effective 
potential V e ff created by the electronic and nuclear charges is given in atomic units as 

y«M-f & + P) 

Zi is the charge of the 7th nucleus located at R. The electron charge density is constructed from 
the Kohn-Sham eigenfunctions tpi as 

oo 

n(r)=^/ i |^(r)| 2 (4) 

i=l 

with fi being the occupation number of the eigenstate 

Many functionals exist for the exchange and correlation energy, the simplest and widely used 
being the local spin density approximation (LSDA) which exists in several forms 000. 
Improvements on the locality are found by using the generalised gradient approximation (GGA) 
(e.g. Perdew et al. 00). 

The index i in the basis set {^(r)} is easily understood in the language of the one-particle 
Bloch states \vk) which are expanded into atomic-like waves inside of spheres centred on atoms 
(muffin-tin spheres) matched to plane waves in the interstitial region. The band state inside of 
sphere s is written in the Slater's augmented plane-wave method as 

oo I 

Ck(r)=£ E DUvk) Ul (e^r)Yl(r) (5) 

1=0 m=-l 

where f = r/r and the radial function ^(e^k, r) is the solution of the radial Kohn-Sham equation 
in the spherically averaged crystal potential. The expansion coefficients D i s m (t , k) are determined 
from the boundary conditions at the sphere's surface. The normalisation condition on iptk( r ) ma Y 
be written as 

oo I 

]T £ I Df m (vk) | 2 = 1, (6) 
1=0 m=-i 

which is a sum over all the local partial charges associated with the band state | z^k). 
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B. Angular momentum projected local density of states and electron energy loss 

spectroscopy calculations 

In the above formulation the angular momentum projected density of states is written as 

xi m (E)^Y,\ D tm(^)\ 2 HE~e^). (7) 

The vector analogy of the p orbitals permits the local p density of states to be projected onto a 
cartesian coordinate system as 

xl, Pv , Pz (E) = ]T |£>^>k)| 2 <5 (E - e „ k ) , (8) 

with D% = -j= [Dl +1 - Dl_ x ],Dl = [23f i+1 + Df,_J and D° z = Df >0 . 

We now introduce the expression central to the theory of electron scattering by the atoms of 
a crystal. The double differential scattering cross-section for the excitation of an atom by fast 
electrons is given in the first Born approximation by |4l| 

where ao is the Bohr radius, 7 = (1 — is the relativistic factor, ko and k the fast electron 

wave vectors before and after interaction, respectively, and Q = ko k is the momentum transfer. 
The dynamic form factor (DFF) S(Q, E) is defined as 

S(Q, E)=J2\ (¥ Q ' T \f) \ 2 5 (E + Ek- E f ) (10) 



for excitations from states \i) with eigenvalue Ei to state |/} with eigenvalue Ef. In terms of the 
coefficients of the augmented plane wave expansion Eq.©, Nclhiebel et al. [4JJ| explicitly wrote 
down the DFF for scattering from an initial state \nlm) to the projection of the final states onto 
the angular momentum state \l'm' } as 



S(Q,E) = jJ^^^M-ir 

I'm' L' M' Xfi A'p/ 



x(2l + 1)V(2A + 1) (2A' + 1) (21' + 1) (2L' + 1) 

XK A (Q)*KY(Q) Ux(Q))ne'W (h'iQ))„e'lL' 



X 



I A l'\ ( I A' L' 
oj 1^0 

I X V \f I X' L' 
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tn 

xJ2Dt>m'(.v*-)Dt, M ,(vkyS (e' - e „ k ) (11) 
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with e' = E + e n i. The radial integrals {jx{Q)) ne , w are defined as 

(h(Q)) ne m>= dRR 2 u nl (R)j x (QR) Ul ,(e',R) (12) 
Jo 

where u n i (R) is the radial part of the core electron's wave function, R t is the radius of the muffin- 
tin sphere and j\(x) is the spherical Bessel function of order A. In the small Q limit the first two 3j 
symbols q boil down to the dipole selection rule I' = l± 1, L' = l± 1. For finite Q. Eq. l|ll|l 
incorporates the orientation resolved character of the ELNES and it can therefore be of interest in 
probing the level of anisotropy present in a crystal and its dependence on the microscope setting. 
The cross terms in the DFF are those corresponding to the coupling of different final angular 
momenta {V / V) resulting from the coherent composition of the band state | i/k). These terms 
become important for large scattering angles. 

For nanotubes, anisotropy effects on ELNES will be studied with respect to the angle 7 between 
the electron beam and the normal to the walls of the tube. In all the orientation resolved ELNES 
calculations, the electron beam, the tube axis and the normal to the walls of the tube will all lie in 
the same plane. The microscope setting will be characterised by the beam convergence semi-angle 
a and the aperture collection semi-angle (3. The dynamic form factor is calculated by integrating 
over all possible orientations of the momentum transfer within the solid angles subtended by the 
convergence and collection semi-angles. Integrated ELNES (relevant for polycrystalline samples) 
will be obtained by averaging over all possible electron beam orientations. 

C. Relevant coordinate system of the nanotube 

On Figure ^ we show the relevant coordinate system of the cylindrical nanotube. The decom- 
position of the local p density of states (pDOS) of Eq.JSJ) for an atom sitting at the point N will 
be allocated as radial (for the r-pDOS), perimetral (for the i^-pDOS) and axial (for z-pDOS). All 
the coordinates of the atoms in the unit cell of a nanotube can be related by symmetry to those 
of a particular atom. Thus, the local electronic properties of all the atoms of the homogeneous 
defectless nanotube should be identical and the total DOS should be a scalar multiple of the local 
DOS of an atom. The decomposition of the local pDOS of an atom into radial, perimetral and axial 
components should therefore reflect that of the total electronic density of state of the nanotube. 
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Figure 1: Relevant coordinate system of the cylindrical nanotube. We shall denote the z direction as axial, 
the azimuthal (ip) as perimetral and the normal to the walls of the tube as radial (r). 



Curvature effects of the tubes on the properties can be studied with respect to the pyramidali- 
sation angle fi. To define the pyramidalisation angle, we consider the zigzag tube (n,0) whose arc 
is drawn on Figure El This arc is depicted by the arcs IP J, EDF and GCH. OZ is along the tube 
axis. On the figure we show a carbon atom D bonded to three other atoms A, B and C. For the 
flat graphene sheet A, B, C and D all lie in the same plane but as curvature is introduced as a 
result of tube formation atom D becomes lifted out of the plane bearing A, B and C. In order to 
quantify the curvature of the nanotube on the local structure of a carbon atom, we may define 
the pyramidalisation angle /i for the 3-coordinated atom D as the angle between the normal to 
the plane bearing the three nearest neighbours (atoms A, B and C) and any of the 3 carbon bond 
vectors (bonded to D, either bond AD, BD or CD). The C-C bond length is fixed at b=1.42A and 

1 /2 

from Figure it follows that z = PD = (b 2 — 2R 2 (1 — cos8 n )) where R is the radius of the tube 
which for an ideal (n,m) tube is R = (n 2 + m 2 + nm) 1 ' 2 . Geometrical considerations lead to 



D. The pyramidalisation angle /.t of zigzag and armchair tubes 



M(n,o) — 7r /2 = arc sin 



(1 - cos(0„)) 



(13) 
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where x n — (40^/3 — 2 (1 — cos(0„))) and 9 n — tr/n. Similar considerations for the armchair 
tubes (n,n) result in 

M(n,n) - k/ 2 = \ arc sin 

To first order in b/R we get /j, ~ 7r/2 + b/(AR) for both zigzag and armchair tubes which is the 
result obtained by Dumitrica et al. 42] . It is important to remark that if we define the tt orbital 
axis as one making equal angles with all three a bonds of a 3-coordinated atom then this angle 
will be equal to the pyramidalisation angle for the case of a nanotube. The deviation of this angle 
from 90° induces a new hybrid between the s and the p z (z in the language of the planar trigonal 
coordination) orbitals (s m p z ). The proportion of s states r = m/(m + 1) in the 7r-orbital has 
been expressed in the 7r-orbital axis vector (POAV1) approximation |43l \A& |45| as r — 2cot 2 /x 
from which the new bonding state is sp 2+ with 2 + 5 = 2(1 — r)/(l — 2r). In other words, a 
fraction 2r/(l — 2r) of the p z orbital participates in the cr-bond. For the Ceo molecule 5 ~ 0.278 
meaning an s fraction of r ~ 0.11 in the it orbital of Cgo- For tetrahedral coordination, fi = 109.47° 
corresponding to a sp 3 hybridisation. 




Figure 2: Geometry for determining the pyramidalisation angle of a zigzag (n,0) tube. OZ is along the 
tube axis and OX and OY are radial to the tube. IP J, EDF and GCH are arcs of the tube. 
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III. COMPUTATIONAL DETAILS 

We study the electronic properties of a series of carbon nanotubes and some curved carbon sys- 
tems using the ab initio all-electron Full-Potential-Linearised- Augmented-Plane- Wave (FPLAPW) 
package WIEN2k 46]. The unit cell is partitioned into two regions; regions of nonoverlapping 
spheres centred on the atoms and an interstitial region. The basis functions inside the muffin-tin 
spheres are linear combinations of spherical harmonics Yj^ (f ) and radial functions ui (r) . These are 
augmented with a plane wave basis in the interstitial region (with the coefficients being functions 
of the reciprocal lattice vectors and the wave vectors of the first Brillouin zone). The maximum 
Z-value of 10 is adopted inside the atomic spheres. The charge density is Fourier expanded up to 
G m ax — 14. The Brillouin-zonc integration is done using the Blochl 47] improved tetrahedron 
method. The exchange and correlation energy is treated using the local spin density approxima- 
tion 37] which is a fit to the Green's- function Monte Carlo calculations of Ceperley and Alder Q]. 
Core states are treated fully relativistically. 

The choice of the muffin-tin radius and the plane wave energy cut-off are governed by the 
parameter RKM = Rmt x K max (where Rmt is the smallest muffin-tin radius in the crystal). 
For these carbon systems the Is state is the core state while 2s and 2p are the valence states. 
Muffin-tin radii of 1.3 a.u. are used. The total number of k-points in the full Brillouin zone was 
fixed at 1000 (between 60 and 80 in the irreducible wedge, depending on the dimensions of the 
unit cell) for the self consistency calculations without the core-hole effect and 300 (between 100 
and 150 in the irreducible wedge - due to lower symmetry) for core-hole calculations in larger 
unit cells. Due to the supercell nature of each unit cell of the nanotube, these values should be 
representative of infinite samplings of the Brillouin zones. The tubes (and also the nonoarcs) are 
arrayed in orthorhombic unit cells. The ideal graphene folded cylinders are adopted for all the 
tubes studied. For the (6,6) tube the difference in DOS and ELNES for RKM values of 5, 6 and 7 
were insignificant. So we adopted an RKM value of 5.5 for all other tubes. 

When a core electron is excited into the unoccupied level a core-hole is left behind which interacts 
with the valence and conduction electrons and modifies the crystal electronic structure if the core- 
hole is not well screened by the valence and conduction electrons. In metals the screening is effective 
meanwhile in semiconductors there is less screening and the core-hole interaction becomes very 
important. This core-hole interaction is a dynamical process and should be treated dynamically 
using the Green's- function approach Two ad hoc techniques are widely used to account 

for this core-state effect. These include the Z+l method or the equivalent core approximation 
whereby the excited core seen by the valence and conduction electrons is equivalent to that of 
an atom with atomic number Z+l. For the case of carbon systems the core atom is replaced 
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by a nitrogen atom. The second method, which we term the core-excited method, consists of 
introducing a core-hole in the core state by removing an electron from it and inserting it as 
a uniform background charge. In the so-called "sudden" approximation no core-hole effect is 
considered and ELNES or any excitation process is calculated from the ground state electronic 
structure. 

Ea. [fl"T|l has been implemented into the WIEN2k package Q- ELNES was calculated on the 
lowest symmetry- adapted k-mesh of each system where all symmetry operations are eliminated 
except the identity. By doing so cross terms of the dynamic form factor are well reproduced. The 
core-hole effect is introduced in all the tubes studied via the core-excited method. For the armchair 
tubes (n,n) the unit cell was doubled along the tube axis in order to increase the distance between 
core atoms. This means a minimum distance between core atoms of 4.9lA and 8n atoms per unit 
cell. Since the length of the smallest unit cell along tube axis was 4. 26 A for the zigzag tube (n,0), 
these smallest unit cells with 4n atoms were used for all such tubes with one of the 4n atoms made 
the core atom. The minimum wall-to-wall distance was set at 4A to reduce intertube interactions. 

IV. RESULTS AND DISCUSSIONS 

The various tubes studied together with their diameters and pyramidalisation angles are shown 
on Table |U In order to study the curvature effect on the electronic properties of nanotubes, we 
also considered the case of the (4,0) tube which has not been shown to exist in nature. The (5,0) 
tube has recently been synthesised in a zeolite cage [3]) . 



Tube 


Diameter (A) 


M°) 


Tube 


Diameter (A) 




Tube 


Diameter (A) 




(4,0) 


3.1 


101.9 


(7,0) 


5.5 


97.2 


(4,4) 


5.4 


97.4 


(5,0) 


3.9 


99.8 


(8,0) 


6.3 


96.2 


(5,5) 


6.8 


96.0 


(6,0) 


4.7 


98.3 


(10,0) 


7.8 


95.1 


(6,6) 


8.1 


95.0 



Table I: The different tubes studied; their diameters and pyramidalisation angles. 



A. Total and local density of states 

In the zone-folding tight-binding approach, any nanotube (n,m) is metallic if | n — m | is divisible 
by 3 26]. From this it follows that all armchair tubes are metallic. Figure shows the ground 
state total DOS of all the studied tubes. 
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Remark the feature at the Fermi level for the strongest curved tube (4,0). This feature is also 
present for all other small zigzag tubes but to lesser intensity as the curvature decreases. It vanishes 
for the (10,0) tube. Such a remarkable curvature dependent feature is absent for all armchair tubes 
whose DOS remains flat within ±1 eV around the Fermi level except for a depression of the DOS 
at the Fermi level for the small (4,4) tube. Such a depression was also reported by Delaney et 
al. [3| on a DFT calculation on the (10,10) tube and this was due to inter-tube interactions. The 
metallicity condition according to the zone-folding tight-binding theory is based on the 7r and 7r* 
bands. The notion of tt and a becomes ill-defined for these curved systems but we may speak of 
curvature induced tt-cj (n and a of the planar graphene) hybridisation to form new hybrids. A 
careful study of the local density of states (for a particular atom) in all three directions (axial, radial 
and perimetral to the tube) reveals some interesting aspects: in that while the finite LDOS at the 
Fermi level of all the small tubes are dominantly radial to the tube (see Figure 0J, non negligible 
axial and perimetral LDOS are present around the Fermi level for the small zigzag tubes. This 
is a signature of ir — a hybridisation. This hybridisation is weak for tubes whose diameters are 
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greater than 5. 5 A irrespective of chirality. This tt — a hybridisation may be quantified in the 
7r-orbital axis vector (POAV1) approximation [43I [ill I45I] described above. In this approximation, 
as much as 22% p z electrons participate in the a bonding for the (4,0) tube while less than 4% 
participates for the (7,0) and (4,4) tubes of similar diameters ^5.5A. It is important to remark that 
while the curvature-induced hybridisation drastically modifies the electronics of small diameter 
zigzag tubes its effect on arm chair tubes is weaker. Blase et al. I&3 studied the hybridisation 
effects on the metallicity of small radius zig-zag NTs and concluded that the curvature induced 
7r — <t hybridisation effects change the energy and character of lowest lying conduction band states 
with important consequences on the metallicity and transport properties of the tubes. Reich and 
Thomsen |28j also pointed out that, compared to the zigzag tubes, the electronic band structure 
is less affected in the armchair and chiral tubes by the curvature induced tt — a hybridisation. 





(8,0) 



(10,0) 



-4 4 
E-E F (eV) 




Figure 4: Axial (full line), radial (dashed line) and perimetral (dotted line) local DOS for nine nanotubes. 
The vertical line denotes the Fermi level. 
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B. Energy loss near edge spectroscopy (ELNES) 

1. Effects of core-hole interaction 

For the (4,0), (5,0), (7,0), (10,0), (5,5) and (6,6) tubes we performed core-hole-included DFT 

calculations, using the excited-core approach, by removing an electron from the Is state and 

inserting it as a uniform background charge. For the (5,0) tube we also introduced the core-hole 

I I 

effect through the alternative Z+l method [5jJ which essentially gave the same result. These 
core-hole techniques have also been applied on graphite and are found to yield the same results 
which are in good agreement with the ELNES measurements |54| up to about 15 eV beyond the 
threshold. 

For the zigzag tubes the minimum distance between core atoms was set at 4.26 A while for 
the armchair tubes it was fixed at 4.91 A. Due to this difference in the inter-core-atom distances, 
we expect stronger intratube interactions between the core atoms of the zigzag tubes than in 
the armchair tubes though this interaction will not strongly influence the C K-edges. Intertube 
core-atoms interactions can be ignored as the least distance between such atoms is more than 7A. 
On Figure El we compare the effects of the inclusion of core-hole on the C K-edge of six tubes. 
To account for instrumental and life-time broadening we convoluted the calculated spectra with 
Gaussians with full width at half maximum (FWHM) of 0.5 eV. The microscope convergence and 
collection semi-angles (a and 0) were set at the magic values (1.87, 3.01 mrad). Orientation 
averaged spectra are shown though at this microscope setting no orientation effect should be seen. 

The immediate effect of the core-hole interaction is the shifting of the 7r* edge to lower energies 
and the redistribution of a* intensity to the lower energy part. A signature of tt — a hybridisation 
is the feature at 5 eV which decreases in intensity as the tube diameter increases from the (4,0) 
tube through the (7,0) tube and is evidenced in the armchair tubes as the shoulders of the lower 
energy part of the a peaks. Comparison with measurements on SWNT [sol I55I reveals that the 
core-hole must be considered if ELNES measurements are to be correctly explained theoretically. 

Knupfer et al. 13 detected a new feature in a single-wall carbon nanotube Is EELS spectrum 
about 2 eV above the ir* peak (i.e. about 5 eV above the Fermi level), which was absent in some 
samples jiSfij . The origin of this feature remained unexplained. On Figure we remark the tt — a 
hybridisation features 2-3 eV above the main it* peak of small tubes. This may suggest that in 
Knupfcr's sample there could have been tubes as small as 4-10 A in diameter. 
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Figure 5: Comparison between the sudden approximation (full lines) and core-hole (dashed line) calculated 
ELNES of (4,0), (5,0), (7,0), (10,0), (5,5) and (6,6) tubes and the curvature effect on the ELNES of these 
tubes. The Fermi level is fixed as the energy origin. 

ELNES spectra were also calculated in the so-called "sudden" approximation for all the nine 
tubes and compared with that of flat graphcnc in Figure El Although it is very clear by now 
that ELNES calculations on nanotubes that do not include the core-hole effects cannot be used to 
explain experimental EELS measurements, such calculations should lead to insights in the ground 
state momentum projected local density of states and the effects of curvature or local deviations 
from planar configurations on the transition matrix elements. As Figure shows, the direct effect 
of curvature is the splitting of the tt* peak into the Ai and A2 (Bi and B2) features for the zigzag 
(armchair) tubes and the splitting of the a* peak into the A3 and A4 (B3 and B4) features for 
the zigzag (armchair) tubes. Suenaga et al. [^J performed the ELNES calculation of (5,0) and 
(10,0) tubes (in the sudden approximation) and reported the splitting of the 7r* peak as curvature 
increases. They also pointed out that as the curvature increases the a* peak decreases. We remark 
also that as the curvature increases the predominantly 7r* peak Ai diminishes while the intensity of 
the hybrid feature A2 grows. The weak overlap between B2 and B3 and the strong overlap between 
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A2 and A3 for the strongest curved armchair tube (4,4) and strongest curved zigzag tube (4,0), 
respectively, demonstrate that B2 is purely of n* character, B3 purely of a* character while A2 and 
A3 are of mixed n* and a* characters. This point emphasises the fact that the curvature induced 
7r — a hybridisation effect is stronger in zigzag tubes than in armchair tubes. From Figure |5] we 
also notice that the core-hole leads to a clear distinction between the 7r* and the a* edge onsets. 




10 20 30 10 20 30 
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Figure 6: Effects of curvature on the ELNES of small armchair and zigzag tubes. No core-hole calculation 
is implemented. Each spectrum is shifted upward for proper visualisation. The vertical lines denote the 
Fermi level. 

2. Energy dependence of the transition matrix elements 

In order to investigate the relationship between angular momentum projected local DOS and the 
calculated ELNES and the effect of core-hole on the transition matrix elements we plot the radial 
(p r ), perimetral (p v ) and axial (p z ) p density of states of an atom together with the unbroadened 
ELNES for the (5,5) tube on Figure [7| The orientation averaged ELNES spectra are calculated 
for microscope convergence and collection semi-angles of 1.87 and 3.01 mrad, respectively. 
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Figure 7: The perimetral LDOS (p^,), the radial LDOS (p r ) and axial LDOS (p z ) of a carbon atom of 
the (5,5) nanotube and the ELNES of the same atom. The upper panel shows the non core-hole-included 
results while the lower panel shows the effect of core-hole. 

The calculated ELNES features are excellently reproduced by the projected DOS for both non 
core-hole and core-hole calculations. The spikes in the DOS are also beautifully reproduced by 
both ELNES spectra. Figure [S] now shows the case of a (10,0) tube. Remark that the agreement 
between LDOS and ELNES is excellent for the momentum resolved ELNES. This strong agreement 
between ELNES features and LDOS shows the weak energy dependence of the transition matrix 
elements though these will be shown to depend strongly on the electron beam orientation and the 
microscope settings. It also shows that the dipole selection rule governs the ELNES spectrum for 
nanotubes because all the ELNES features are those of the pLDOS. Knupfer et al. [20] reported 
the ELNES spectrum of a SWNT which was free of singularities and concluded that the core- 
hole effect was responsible for this. Though their measurements were done using a high energy 
resolution microscope (0.115 eV), which should be able to detect the spikes, these were conspicuous 
by their absence. Thermal broadening could have been responsible for the washing out of these 
singularities. 
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Figure 8: Local density of states (LDOS) of the carbon atom of the (10,0) nanotube and the ELNES of 
the same atom with core-hole included: momentum transfer orientation resolved calculations is shown. 



3. Orientation resolved ELNES 



On Figure we plot the total ELNES of (6,6) and (10,0) tubes and graphite for two incident 
electron beam orientations. The electron beam makes an angle 7 with the normal to the graphite 
plane or the walls of the nanotubes. When 7 = 0° momentum transfer is parallel to the walls 
of the tube and scattering is mostly due to the in-plane bonds. For 7 = 90° momentum transfer 
is normal to the walls of the tube and scattering probes off-plane bondings. Both the (6,6) and 
(10,0) tubes spectra are for core-hole calculations and are compared with core-hole calculated 
graphite. Orientation resolved ELNES calculations for the (6,6) and (10,0) tubes reveal that 
the magic convergence and collection couple (a*,/3*)=(1.87,3.01mrad) found for graphite remains 
magic for these curved carbon systems (and possibly for all other carbon materials) as Figure [5] 
shows. Anisotropy is seen for j3 ^ (3* when a* is fixed at 1.87 mrad. For (3 < (3* Is— > a* transition 
is strong for the electron beam perpendicular to the graphene plane (7 = 90°). This tendency 
reverses for /?>/?*, all in similar trend as in graphite. Stephan and co-workers 56] measured the 
ELNES for the electron beam normal (7 = 0°) and perimetral (7 = 90°) to a multi-wall nanotube 
and found that the anisotropy was weak compared to that of graphite. A look at the topmost 
two panels of Figure |5] confirms this finding. Following the discussion of the previous section, we 
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can conclude that this anisotropy is entirely governed by the transition matrix elements and that 
the latter depend strongly on the electron beam characteristics and the microscope settings in a 
manner well elaborated in refs. 57. 581. 




Figure 9: Orientation resolved ELNES of (6,6) and (10,0) tubes (bold line) and comparison with graphite 
(thin line). The spectra for 7 = 0° (dashed line) and 7 = 90° (full line) coincide exactly over all energies 
for /3 = 3.01 mrad but no such coincidence is seen for /3 = 6.28 mrad and /3 — 0.90 mrad. The vertical line 
indicates the Fermi level. 

We performed a (l',m') decomposition (which consists of the various terms of Ea. 1)11(1 when 
I' and m' are fixed) of the ELNES in a fixed reference frame. In this decomposition, the Is— » 2p 
ELNES has three components: (1,0), (1,-1) and (1,+1) components. In the case of a planar system 
where the z-axis is perpendicular to the plane (for example, graphite) the (1,0) component will 
correspond to the Is— ► it* transition while (1,-1)+(1,+1) will correspond to the Is— » a* transition. 
In the case of the (6,6) and (10,0) nanotubes, where the z-axis was along the tube axis, the (1,0) 
component becomes a transition into the axial a* orbitals. Unfortunately the (1,-1) and (1,+1) 
components become inseparable mixtures of 7r* and the perimetral a* components for these cases 
(see Figure ITU)) . It is important to point out that the ELNES spectra on Figure ITUI are entirely 
given by the three p ELNES-components for all the three collection angles. This illustrates the 
fact that even for relatively large collection angles, the dipole sellection rule is still valid. In the 
case of the (5,0) tube we fixed the z-axis normal to the walls of the tube such that the (1,0) 
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component can be obtained separately. On Figure El we show the total ELNES, the ir* and 
the a* components of the core-hole calculated ELNES of the (5,0) tube. Remark that, as a 
signature of 7r — a hybridisation, the tail of the a* component overlaps with the tt* spectrum. 
This is in agreement with the observation (see Figure 0} that around the Fermi level there is a 
non negligible DOS resulting from perimetral and axial a orbitals for the (5,0) tube. According 
to the POAV1 approximation (43, |4J, |45j , 14% of the p z orbitals should participate in the hybrid 
orbital formation for the (5,0) tube. This agrees with the 15% obtained when two Gaussian centred 
functions, one (Pi) describing the overlapping lower energy tail of the a* edge and the other (P2) 
for the main a* peak, are used to fit the a* edge as shown on Figure ITTI The 15% was given by 
the ratio Area(Pi)/Area(Pi+P2). Orientation resolved ELNES measurements should be able to 
reveal this rehybridisation in such small tubes. This technique of (l',m') decomposition has recently 
been applied to hexagonal BN 59| and graphite [541 ], both of which have planar symmetry. We 
believe that this w* [a* decomposition in planar systems is a viable alternative to the existing 
techniques for sp 2 /sp 3 quantification of carbon materials based on the isolation of 

the 7r* / a* excitation cross-sections. 




Figure 10: Spectral (l',m') decomposition of the ELNES spectra of the (6,6) and (10,0) tubes. The total 
spectrum (full line) for 7 = 0° are decomposed into the (1,0) or axial component (dashed line) and the 
(1,-1) + (1,+1) or radial+perimetral terms (dotted line). 
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Figure 11: Spectral (l',m') decomposition of the ELNES spectrum of the (5,0) tube. The total spectrum 
(full line) for 7 = 0° are decomposed into the n* and the a* components. Remark that as a signature of 
tv — a hybridisation, the tail of the a* component overlaps with the n* spectrum. A Gaussian broadening 
of 1.8 eV FWHM is used. 



C. Investigation of the local nature of ELNES 



It is usually postulated that ELNES probes the local properties of materials. In particular, 
spatial and momentum resolved ELNES are often used as fingerprints of bonding states and local 
disorder or to explore the presence of grain-boundaries in materials. Due to the strong curvature 
applied to graphene to form a nanotube, carbon can undergo a wide range of bonding states from 
sp 2 for flat graphene to sp 2 ' 28 for tubes with similar diameters as the Ceo molecule. This may 
imply a curvature induced re-activity of nanotubes different from that of graphite 

In this section we consider arcs of tubes which are repeated in space periodically as bonded or 
wavy arcs in order to investigate the local nature of the energy loss near edge structures. We 
calculated the DOS and ELNES of these arcs and compared them with those of the full tubes. 
Stcphan et al. simulated the curvature effects on the unoccupied density of states of NT's 
by considering a corrugated graphene sheet in which the corrugation was defined by the chair 
configuration of the carbon hexagons with a corrugation angle of 8° (an atom is lifted out of 
the plane such that the lifted bond vectors make an angle of 8° with the plane) to model a 10 
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A tube. Such a model which considers only the corrugation angle and ignores the chiral angle 
cannot be expected to be representative (electronically) of tubes as small as 10 A in diameter. 
Our bonded/ wavy arcs are aimed at addressing this problem. 

We define arcs of a given nanotube, e.g. an arc of an armchair tube (n,n) as a section of the tube 
that is cut out and is repeated in space (because of the use of the periodic boundary conditions) 
either by bonds (of the order of C-C bond length - see Figure IT^l or as the wavy objects shown on 
Figure El This latter model diminishes the unphysical hybridisations that are introduced by the 
boat configurations (for armchair arcs) and the 8 atoms ring configurations (for zigzag arcs) in the 
bonded arcs. The size of each arc is chosen large enough so that at least one of its carbon atoms 
is situated at a minimum distance of two C— C bond-lengths from the point of contact of two arcs. 
This is to ensure a reduced interaction between the defectuous points and the carbon atom whose 
local properties should be comparable to those of an atom in a defect-free tube of similar curvature 
and chirality. 

On Figures El an d El we show the total density of states and the ELNES spectra of the C 
K-edges of the (6,6) and (10,0) tubes, respectively. For the ELNES spectra on these figures no 
core-hole calculations are used. The DOS and the carbon Is excitation edges are also shown for 
the nanoarcs. The wavy and the bonded arcs reproduce qualitatively well (up to 12 eV beyond the 
edge onset) the ELNES of the (6,6) and the (10,0) tubes. These good agreements of the calculated 
ELNES of curved arcs with those of tubes of similar radius of curvature and chirality confirm the 
fact that ELNES (up to 10-15 eV above threshold) probes local structures. The total density of 
states, on the other hand, is a property that depends on the underlying periodicity. It should be 
strongly affected by periodically repeated defects like the boats and 8 atoms rings present in the 
bonded arcs. While the LDOS describes the local electronic structures, the total DOS is responsible 
for the observed macroscopic properties like metallicity. 



Figure 12: Armchair and zigzag bonded arcs. 




Figure 13: Zigzag and armchair wavy arcs. 
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Figure 14: The total DOS and ELNES of the (6,6) tube and its two model arcs. 




Figure 15: The total DOS and ELNES of the (10,0) tube and its two model arcs. 
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V. CONCLUSIONS 

We have performed first principles calculations on a series of armchair and zigzag nanotubes 
of diameters less than fnm using the all-electron Full-Potential(-Linearised)-Augmented-Plane- 
Wave method. We have investigated in detail the effects of curvature, core-hole interaction and 
orientation on the ELNES of carbon nanotubes. The inclusion of the core-hole effect is very 
important to reproduce the measured ELNES spectra of carbon nanotubes jsoi l.~>(| . A curvature- 
induced 7r — (j hybridisation is pointed out and shown to have a stronger effect on zigzag tubes 
than on armchair tubes. Not only does ELNES calculated in the sudden approximation [2jj fail 
to explain the experimental measurements, an analysis of the tt — a hybridisation based on it may 
be misleading especially for small diameter zigzag tubes as they do not show a clear separation 
between the tt* and a* features. We have shown that the ELNES of carbon nanotubes obeys the 
dipole selection rule and that the local DOS of the excited atom and the ELNES spectrum may 
show spikes due to the Van Hove singularity, which are characteristic of the quasi one-dimensional 
nature of the nanotubes. The curvature effect is translated into the electron energy loss spectra by 
the splitting of the tt* and the a* edges. A reduction in the anisotropy, as compared to graphite, 
is also recorded for these highly curved carbon systems in line with the findings of Stephan et 
al. although this anisotropy depends on the microscope settings. It is found that the magic 
convergence and collection semi-angles (1.87, 3.01 mrad) found for graphite gives an electron beam 
orientation independent ELNES for nanotubes. Therefore, nanotubes can be used as an alternative 
to graphite as an anisotropic material to determine the magic collection and convergence semi- 
angles for carbon materials. The advantage of nanotubes over graphite is the ability to explore 
a large range of focused electron beam orientations. The (l',m') ELNES decomposition has been 
shown to be very valuable to unambiguously quantify the tt — a hybridisation in small nanotubes. 
This decomposition may be a promising technique to study the type of bonding in carbon materials. 
We have performed calculations on model arcs of tubes and confirmed that ELNES probes the local 
bonding structures of materials. 
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